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ABSTRACT 


A proof-of-concept experimental validation of a proposed idea for 
a prototype electro-optical phase encoded position transducer was 
conducted. The intensity of two IR LED beams were agduiated 
sinusoidally at frequency o, with a 90 degree temporal phase 
dizference’ They were further modulated by a pair of sinusoidal 
Jepeteat encoding masks with 90 degree spatial phase difference. 
The mask pair was mounted on a mechanical stage and translated 
perpendicular to the beam axis. The sum of the two signals 
produced by this electro-optical configuration constituted another 
sinusoid at temporal frequency w, whose phase was proportional to 
the mask position. Although small deviations from ideal behavior 
were observed, the validity of the technique has been established 


incontestably. 
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I. INTRODUCTION 


A. MOTIVATION 

This thesis project is part of a larger research and 
development program to design and build a next generation 
middle infrared detection and imaging system by the author's 
advisor, D. S. Davis. Part of Davis' proposed technique is 
that a state-of-the-art imaging spectroscopic system can be 
dsyaleved using orthogonal Walsh functions to _ encode 
multiplexed images, thereby achieving extreme sensitivity 
[Ref. 1]. The validity of this technique was proven by R. H. 
McKenzie, Capt., USMC, in an earlier thesis [Ref. 2], 
and further investigated in a thesis by B. J. Musselman, LT., 
USCG, [Ref. 3]. 

This project stems from Davis' interest in spectroscopic 
research. He has previously designed a cryogenically cooled 
infrared multiplexing Fourier transform epectromsters which is 
flown on NASA's Kuiper Airborne Observatory (KAO). Although 
the new instrument was initially proposed as a replacement for 
the current instrument to use on NASA's next generation 
airborne observatory SOFIA (Stratospheric Observatory For 
Infrared Astronomy), the technology should also be useful in 
other applications such as passive target surveillance 


detection and identification. 





B. THESIS GOAL 

Davis' technique requires the precise mechanical 
positioning of a pair of focal plane encoding masks, by means 
of an electro-optical servo system. Experience with previous 
generations of multiplexing instruments has determined that 
the best servo system ‘seeeovnnnea “ie achieved with servo 
feedback signals that. are phase-sensitive functions of 
loop error, rather than simple magnitude functions 
[Ref. 4, Ref. 5]. “the goat of thie thesis is to develop and 
to implement a proof-of-concept. validation for such an 
electro-optical phase encoded position monitor, as it will be 
employed in the overall multiplexed imaging system. 

The proposed system will be sensitive to the relative 
phases of a pair of amplitude-modulated IR beams. The 
intensity of the beams will be sinusoids with a 90° temporal 
phase difference, i.e., in quadrature, and modulated at 
frequency o. They will be further modulated by a pair of 
sinusoidal optical encoding masks with 90° spatial phase 
difference. Let x represent the position of the masks. Then 
the optical configuration will produce two signals, one 
proportional to sin(®#t)sin(x) and the other to cos(®t)cos(x). 
Their sum constitutes another sinusoidal signal at frequency 
wo, whose phase is just the mask position relative to some 


origin, 


cos (@t) cos (x) tsin(®@®t) sin (x) «cos (®@t-x). (1) 








Careful monitoring of the temporal phase difference between 
this signal (at frequency w) and the original phase reference 
signal (also at frequency »), should yield a phase difference 
that is linearly proportional to mask position. This system 
would allow for extremely accurate measurements of the mask 
position due to the precision with which phases can be 
measured. The position measurements would also have the 
advantage that they would not have the susceptibility to beam 
amplitude noise, such as that produced by power source 
fluctuations. 

Figure 1 shows a block diagram overview of the project 
design. The electronics section consists of four main 
modules. The light sources module produces the quadrature 
sinusoid signals used for driving the LED emitters. The 
detector circuits convert the optical signals modulated by the 
optical masks to electronic signals. The signal processing 
circuits then sum the two output channels and filter out 
nonessential signal elements. Finally the phase comparator 
measures the relative phases between the reference (driving) 
signal and the summed output signal. The optical section 
consists of four positive lenses to focus the LEDs' output and 
the mask pair that modulates the light signals. Each of these 


elements is discussed in detail in the following chapters. 
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Figure 1 Block Diagram Overview of Proposed Thesis Project 








II. ELECTRONICS 


The electronic portion of the experimental setup consists 


of three parts, listed below and described in detail later in 
this chapter. 

e Two-phase quadrature sinusoidal oscillator circuit. 

e Optical detector circuit. 


e Signal processing circuit. 


A. OSCILLATOR CIRCUIT 
; 1. LED . 

Drivers for the LED light sources were the first 
circuits to be constructed. Initially it was not known what 
modulation frequency would be required to achieve the best 
results, since the frequency bandwidth characteristics of the 
LEDs and detectors were not specified. In order to determine 
the optimum frequency experimentally, it was decided to drive 
a standard IR LED (peak wavelength A=950nm) with a Tektronix 
CFG250 function generator and to measure the LED's frequency 
response. IR light sources were chosen instead of visible to 
allow working with the room lights on ana to avoid 60 cycle 
pick-up from the room lights. It was found that the function 
generator had insufficient current drive capability to power 


the LED. Therefore, the output of the function generator fed 





an emitter follower, which, in turn drove the LED. An emitter 
follower current amplifier, using a type 2N3053 NPN silicon 
transistor (component Q1 in Figure 2), was designed and 
constructed [Ref. 6:pp. 68-71]. 

The supply voltage for the emitter follower is +5V and 
is provided by a Texas Instruments uA7805c positive voltage 
regulator. The light source itself is an Archer SY-IR53L high 
intensity infrared LED. The LED requires a maximum sustained 
forward current of 20mA anda forward voltage of at least 1.3 
vpc. This means that, theoretically, an external emitter 
resistance of at least 65n is needed for the LED to light 
without drawing aeceasive current. . With the resistance used 
(Ry, in Figure 2) the measured voltage drop across the LED was 
the required 1.3 vpc. 

It is necessary to biax the follower's base to ensure 
that a collector current flows during the entire input signal 
swing. The simplest way to do this is with a voltage divider 
[Ref. 6:p. 70]. The 1.3 VDC drop across the LED combined with 
the internal 0.7 VDC drop across the transistor mandates that 
the base voltage be at least 2.0 VDC above ground. The 
negative swing of the source signal must aes be accommodated, 
so the base voltage actually needs to be higher than 2.0 VDC. 
The values of the bias dividing resistors, R,, and R,,, were 
chosen to keep the measured base voltage at 3.4 VDC. T h e 


input resistance R,, and the load resistance R,, are chosen to 























Light Source Schematic Drawing 


Figure 2 





be sufficiently small so that the effective impedance of the 
emitter follower is not established by the base biasing 
resistors. - 
2. OSCILLATOR 
The two-phase oscillator that drives the LEDs in 
quadrature has to operate in the 70kHz to 125kHz frequency 
range. This range was determined by testing the frequency 
response of the LED/detector combination, as explained below 
in section II.B.3. Although there are several ways to produce 
two signals with the desired phases, a two-step method was 
chosen. The first step is to generate quadrature square waves 
using digital circuitry. The second step is to filter out the 
square waves! overtones using a low-pass filter, so that only 
the fundamental Fourier component, at the desired modulation 
frequency, remains. The details are described below. 
a. Timer 
The first component in Figure 2, U1, is the square 
wave janevetor a Texas Instruments TLC555C timer chip, the 
operation of which is straightforward. Initially, when power 
is applied to the chip the output is HIGH, the charge on the 
capacitor is low, and the internal discharge transistor is off 
{Ref. 7:p. 171]. The capacitor then begins to charge 
through the series resistors R, and R, toward the supply 


voltage of +5V. When the capacitor gets to % V,. (or 3.33 V) 


the threshold is triggered, causing the output to go LOW and 








the discharge transistor to turn on [Ref. 6:pp. 286-287]. The 
capacitor is now discharged through resistor R, until the 
capacitor voltage drops to % V, (1.67 V). At this point the 
threshold triggers and returns the chip to its initial state. 
The operation is now cyclic with period T=0.693(R,+2R,)C. For 
the values shown, this gives a calculated frequency of 306kHz, 
and a measured frequency of 312kHz [Ref. 7:p. 47]. The 
nengured frequency is approximately four times the desired 
operating frequency. This is necessary because subsequent 
production of the quadrature signals involves four-fold 
frequency division. This is described in detail in the next 
section. 
b. Frequency Divider and Phase Shifter 

The quadrature phase shift and frequency division 
are accomplished using two National Semiconductor DM74S74AN 
flip-flops (shown in Figure 2). 

The rectangular pulse train A in Figure 3, at four 
times the desired frequency, is input from the TLCS55C 
oscillator to the clock (CLK) input of U2A in Figure 2. 
Because the NOT Q output of U2A is tied to the D input, each 
positive leading edge of pulse train A toggles U2A to its 
alternate state. This produces the 50% duty cycle square wave 
B (Figure 3), at half the input frequency or twice the desired 


LED modulation frequency. The Q output is then fed to the 


clock of U3B. U3B is configured the same as U2A, so that its 











Q output is the square wave C in Figure 3, at the desired 
modulation frequency. The NOT Q output of U2A is connected to 
the clock input of U3A, and therefore U3A responds to the 
negative aetna transition of pulse train B by passing the 
compliment of C to the Q output of U3A. However, the 
frequency of B is twice that of C, so that its negative going 
edges are phase shifted by one quarter cycle. This forces the 


square wave D (Figure 3) to lead pulse train C by 90 degrees. 


Figure 3 Square Wave Pulse Trains 


c. Filter 
Once quadrature square waves have been generated, 
the next stage is to filter out the overtones to produce a 
sine wave. This can be accomplished by passing the signal 
through a low pass active filter with a cutoff frequency, f,, 


between f,,, and 2f,,, where f,,, is the desired LED modulation 


frequency. The Butterworth filter configuration was chosen 
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because it gives the flattest passband possible, combined with. 
a reasonably fast initial falloff and moderate overshoot 
[Ref. 8:p. 72]. The design for the filter was adapted from 
Horowitz and Hill [Ref. 6:pp. 274-275] and consists of ICs U4A 
and U4B in Figure 2. 

By cascading two filter sections together, the 
performance of the filter can be improved. The two sections 
of the filter are not identical but instead are designed so 
that each one contributes by adding its own cutoff frequency 
and damping as a quadratic factor of an nth-order polynomial 
describing the overall filter response [Ref 8:p. 18]. 

The frequency response properties of the filter are 
set by the input resistors R, 5,8,9 and capacitors Cy 56,77 whose 
values are determined by f.=1/2mRC. Resistors Re ip are 
determined from R=(K,-1)R;, where K,=1.152, k,y=2.235 
{Ref. 6:p. 274]. The resistors at the output of the op-amps 
provide a non-inverting voltage amplifier with gain 
KA=1+Re to/Rz 41 where n is the filter stage. The values shown 
in the figure are for a cutoff frequency of £,=100kHz. A 
National Semiconductor LF412 dual op-amp was used so that both 
stages of the filter could be included on one Ic. 

d. Offset Trim 

The output signal from the two-phase oscillator is 

not yet sufficient to drive the LEDs. The characteristics of 


the signal and its DC offset do not match the output signal 
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from the function generator that drove an LED during the 
initial prototyping and testing discussed above in section 
A.1. By measuring the output of the function generator, as 
used during initial tests, it was concluded that the driving 
signal must have a voltage of 3 Ven with a DC offset of 3.5 V. 
This was necessary to ensure operating in the linear range of 
the LED I-V curve. The output from the oscillator is 4v,, with 
a 2 vbc offset. The wietege aay to match the oscillator's 
output characteristics with the emitter follower requirements 
was to insert a non-inverting amplifier with a voltage divider 
at the summing point (U6A) in Piguke 2. The resistors Ris 16 
divide the applied -15 VDC signal by two. IC U6, which is an 
LM412 dual op-amp so that both channels can be accommodated on 
one chip, then amplifies the voltages with a gain of three, as 
determined by Ris and R,.. These two signals are added 
together and the voltage is further divided down to the 
required amplitude and offset by resistors Rig 20° 
| A unity gain amplifier (U5, Figure 2), using the 
LM310 unity gain buffer op-amp immediately upstream from the 
filter, isolates the amplifier from the oscillator. The 
resistor R,, is recommended by the manufacturer. 
e. Completing the Circuit 
After the successful prototyping of the two-phase 
oscillator, a more permanent version of the circuit was 


assembled on soldered circuit boards. The ICs were attached 
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to the circuit boards uSing solderless sockets. Pin-to-pin 
connections were made using above board soldered wire jumpers. 
The circuit boards were then screwed onto a plywood base with 
aluminum standoffs. The power supply leads from the circuit 
boards were connected to terminal blocks attached to the 
plywood base. A photograph of all the completed circuit 


boards is provided in Figure 4. 


B. DETECTOR CIRCUIT 
1. INITIAL DESIGN 
The next phase of the development was the design of a 
matched pair of detector circuits. The main component of each 
circuit is a PIN photodiode. The initial design was simply to 
connect the detector to an inverting transimpedance amplifier 
using a LM301 op-amp. A gain greater than unity was required 
in order to ensure an output signal large enough to use 
downstream. The values of the input resistor R, and the 
feedback resistor R, were chosen so as to provide ample gain. 
2. PROBLEMS/SOLUTIONS 
The photodiode used for the detector has an internal 
impedance Z~10°N when operated in the kilohertz range, 
increasing to Z~108n at lower frequencies. This allows a very 
small photocurrent to produce a large voltage (Ohm's Law). 
Consequently, the circuit is highly susceptible to unwanted 
electrostatic coupling as well as 60Hz pickup (magnetic 


coupling) from its surroundings [Ref 6:pp. 455-457]. 
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Photograph of Completed 


Figure 4 








A shielded box was built to help minimize coupling 
problems. The box borers of an aluminum casing with two 
apertures for the photodiodes and BNC connectors for the 
output. An internal aluminum plate placed immediately above 
the circuit board completes the shielding. 

A shielded box however, does nothing to reduce the 
magnetic coupling. By using shielded twisted pair wire to 
connect the photodiodes to the op-amp, trimming the leads of 
the photodiodes as short as possible, and using twisted wire 
pairs to connect the op-amp output to the BNC connectors, 
pickup and interference are reduced to acceptable levels. 
Replacing the power supply with four 9V batteries, connected 
in series to provide +18V to the circuit, added additional 
protection. 

As a further precaution, a unity gain amplifier (U2, 
Figure 5), using an LM310 op-amp, was added between the 
inverting gain amplifiers and the BNC connectors ‘to isolate 
the high impedance of the detector circuit from the downstream 
components. 

3. TESTING 

Once the circuit was completed, the next step was to 
test the LED/detector circuit to determine the best operating 
frequency. The detector output was connected to a Hewlett- 
Packard Fourier analyzer, while the LED was driven by a 


function generator which was set to frequency sweep mode. 
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Figure 5 Detector Circuit Schematic Drawing 
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Overall LED/detector sensitivity was measured as a function of 
modulation frequency. This was necessary because the internal 
resistance and capacitance of the detector and the AC coupling 
effects of the low-pass filter combine to give an optimum 
frequency range. The results of the Fourier analysis 
established that the best operating frequency for the system 


was in the range between 70kHz and 125kHz. 


C. SIGNAL PROCESSING 
1. SUMMING AMPLIFIER 

The main component in the signal processing section 
(Figure 6) is the summing amplifier (U1) that adds the two 
detector output signals in accordance with Equation (1). The 
input signals at this point must have the same amplitude in 
order to prevent scaling imbalance problems with the 
trigonometric identity (equation (1)). meratere: 
potentiometers are used for the input resistors R, and R, so 
that the gain of each signal can be adjusted as necessary. 

The procedure for matching the signal gains was 
simple. With the input for channel two (the sin(ot)sin(x) 
signal) grounded, the op-amp acted as a simple inverting gain 
amplifier for channel one (cos(wt)cos(x) term) with gain of 
-R,/R,+R,. The amplitude range was measured as R, was adjusted 
through its full range. The procedure was repeated for 


channel two. An output amplitude that was within both ranges 
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Signal Processing Schematic Drawing 


Figure 6 
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was then chosen, and R, and R, were adjusted to provide the 
necessary gain. | 
2. HIGH PASS FILTER 

The output of the above circuit needs to be further 
processed: there is DC term as well some low frequency terms 
that need to be filtered out. This was accomplished with a 
high-pass Butterworth filter. This filter is similar in 
design to the low-pass filter described previously, with the 
only difference being an interchange of resistors sie 
capacitors. For this filter a two-pole design with any cutoff 
frequency below the modulation frequency f,,, will produce an 
acceptable response. The actual cutoff frequency chosen is 
£,=15kHz. 

3. FINAL PROCESSING 

The final two components (U3A and U3B in Figure 6) in 
the signal processing circuit are summation amplifiers used to 
extract cos(ot) and sin(wt) terms from the signal output from 
U2. IC U3A has as its tipute the two reference signals 
driving the LEDs, with two potentiometers as the input 
resistors R, and R,,. The potentiometers are adjusted for the 
correct gain and the voltages are summed. This signal 
cos(wt)+sin(ot) is then added to the output of the high-pass 


filter. The second summation amplifier is set for a gain 


of -1. 








III. ‘OPTICS . 


A. OPTICAL LAYOUT 
1. COMPUTER GENERATED MASK 
a. Mask Design 
Once the electronic circuits were complete, the 
next phase was the dptical design. The first optical problem 
was the design of a mask pair to modulate the two signal 
channels sinusoidally as a function of the mask's position. 
Both channels must be modulated simultaneously and in 
quadrature as the mask is translated. 
| The initial mask design was a sine wave spaced 
vertically 5cm above a cosine wave, with transparent areas 
between the curve and the x-axis and IR-opaque areas "outside" 
the curve. It was decided, however, that this aesign would 
not work because the detector would not recognize when the 
sine wave (or cosine) went negative. For example, as the mask 
was moved from x=0 to x=m, the intensity of the light that 
passed through the mask would increase from zero to a maximum 
value and then to zero again. As the mask moved from x=7 to 
x=2n the intensity would follow the same pattern. In other 
words, the intensity would go from a minimum to a maximum back 


to a minimum in one-half cycle instead of one cycle, because 
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the detector signal was proportional to the inherently 
positive intensity of the detector. 

| A second mask was designed to give a DC offset by 
plotting the curve A(1+sin(wt)), where the constant A is the 
amplitude of the curve. By allowing light to pass through the 
aperture defined by the curve and the x-axis, the intensity 
increases from a mastmam value at x=0 through a minimum value 
and back to the maximum value at x=27. A similar plot for the 
curve A(1+cos(ot)) is drawn 5cem below the sine curve. An 
example of this design is presented in Figure 7. A pair of 
slits, vertically aligned, is then used to reduce the apparent 
area of the mask seen by the detector. 

b. Computer Configuration 
The program to create the optical mask described 

above was written on an IBM Personal System/2 Model 70. 386 
computer. The prototype encoding masks were drawn on a 
Hewlett-Packard model 7550A pen plotter connected to the IBM 
computer via a conventional 9600 baud RS-232 interface. 

The C programming language was selected for the 
program development, based on the author's familiarity with 
that language. An additional contributing factor in the 
choice of language was the availability of the IBM-PC version 


of Turbo C++ 2.0 by Borland International. 
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Figure 7 Optical Mask Design 





c. Producing the Mask 
The copy of the program that produces the mask is 
included in Appendix A. When the program is run, the 
amplitude of the desired curve and the period are specified by 
the user. The program then determines the solutions for the 
equation y=sin(x) and y=cos(x), and passes the solutions to an 
output file called MASK.DAT. This file is in Hewlett-Packard 


Graphics Language (HPGL) format, consisting of two-letter 
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mnemonic codes followed by numerical parameters. The overall 
size of the generated mask is 11x5 cm, as dictated by the 
dimensions of the film holder used to hold the mask in the 
optical setup. Once the mask was produced on paper, it was 
transferred onto a transparency for use during the experiment. 
2. OPTICAL SETUP 

The initial idea for the optical layout was to have 
two positive lenses with the same focal length f=10cm coaxial 
with each mask, LED and detector, as shown in Figure 8. The 
first lens would focus the light from the LED (located at 2f) 
down to a point on the mask (at a distance of 2f) and the 
second lens, also located a distance 2f from the mask, would 
refocus the light onto the detector. A narrow slit would be 
placed as close as possible to the mask to restrict the range 
of x values illuminated by the beams. This idea was rejected 
without testing because the amplitude of the sine wave on the 
mask would have to be extremely small for the focused light to 


cover the whole vertical dimension of the aperture. 


<a 
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LED Ll MASK SLIT L2 DETECTOR 
Figure 8 Initial Optical Layout 
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Figure 9 shows the final version of the layout. 
lens L, is located a distance f from the LED so that the light 
passing through L, is collimated. The light then passes 
through the mask, slit and L, (distances S,, S, and S. are 
arbitrary). Finally, L, focuses the light onto the detector 


located at a distance f. 


<[_ 
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Figure 9 Final Optical Layout 


The optical system was configured with the two channels 
_stacked vertically, with channel one located 5cm above 
channel two. Figure 10 shows a photograph of the final 


optical layout. 
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Iv. RESULTS 


A. INITIAL TESTING 
1. Problems 
Once the experimental setup was complete, the system 


was initially tested by connecting oscilloscope leads to the 


output of the voltage adder (IC U1, Figure 5) and to the 


reference signal. The first of these is the desired phase- 
encoded signal, while the second is the phase reference 
signal. The relative phases of the two signals were then 
visually compared as the position of the mask was changed. 
The results of this test were disappointing. It was 
observed that the pelative phase of the output signal shifted 
to a maximum of ninety degrees, as expected, when the mask was 
translated by one quarter wavelength. However, it then 
shifted back to zero instead of Soneinuens to increase 
linearly. During discussions of possible causes for this 
behavior, the author and his advisor realized that the 
experimental configuration did not account for the DC signal 
terms caused by the fact that light intensities cannot be 
negative. Thus the two channels' intensities are actually 


proportional to 1+sin(wt) and 1+cos(wt) instead of sin(ot) and 
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cos(wt). Consequently, the summed channel signals actually 


generate a voltage of the form 
TeDC+sin (wt) +3in(x) +cos (wt) +cos (x) +cos (wt-x) , ~ (2) 


rather than the one required by Equation (1). 
2. Solution 

The solution to this problem is straightforward. To 
eliminate the unwanted DC aiid low frequency (sin(x) and 
cos(x)) terms in Equation (2), a high pass filter was added 
downstream from the summation amplifier (see section II.C.2.). 
To eliminate the cos(ot) and sin(wt) terms, a second summing 
amplifier was included that would adjust the gain of the LED 
driving signals and add them together before subtracting that 
signal from the high-pass filter output signal, as discussed 


in section II.C.2 and shown in Figure 6. 


B. NEXT ATTEMPT 

With the above changes to the circuit completed, the 
system was ready for the next experimental trial. This time 
a Stanford Research Systems SR510 Lock-In Amplifier was 
included in the setup to monitor the relative phases. The 
output from the circuit wae input into the lock-in amplifier 
as was the reference signal. The initial phase offset between 
the signals was nulled and a lock-in phase shift of ninety 


degrees was added in order to facilitate subsequent phase 
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measurements. The mask was moved in 0.1 inch increments, and 
the resulting phase differences were recorded. 
1. Problem 

The data from this experimental run are listed in 
Appendix B and plotted in Figure 11. The plot shows the 
observed relative phase difference (¢) vs the mask position 
(x). The straight line is a least squares fit of the data,. 
and the curved line is included to show the periodicity of the 
data. This plot shows that the phase difference does increase 
with mask position, as expected, but there is wree a 
significant periodic perturbation from linearity. The period 
of the perturbation is approximately equal to the period of 
the mask. | 

2. Solution 

Possible causes for this periodic anomaly were then 
investigated. The initial belief was that it was due to a 
gain mismatch between the two signal channels. If there was 
indeed a gain mismatch, the summed signals would be 
Acos(wt)cos(x) and Bsin(ot)sin(x), with A¥B, and the trig 
identity (Equation (1)) would not work. Hence, it was decided 
to check for and to minimize any disparity with the gains. 

The first place that the gains were adjusted was at 


the voltage adder IC Ul, Figure 6. The output signal from 
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Figure 11 Phase vs Position for the Initial Data Set 


detector channel one’ was connected to the voltage adder and 
the adder's channel two input was grounded. The amplitude of 
channel one was measured with the oscilloscope and recorded. 
Then channel two was connected to the detector, channel one 
was grounded and the amplitude of channel two was adjusted to 


match that of channel one. 


‘Throughout this thesis, channel one is always the top 
optical system, and is considered to be the cos (wt) cos (x) 
term, while channel two is the sin(wt)sin(x) term. This 
distinction is completely arbitrary. 
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The only other point in the circuit where a gain 
mismatch could possibly be a problem was at the voltage adder 
IC U3A in Figure 6, where the reference signals, sin(wt) and 
cos(wt), are subtracted out (Equation (2)). Channel two in 
the signal processing circuit was initially grounded, and only 
one reference signal was connected to the adder. The output 
of the adder was connected to the oscilloscope. The position 
of the mask was adjusted until the signal measured on the 
scope was at half its maximum value, at which point the mask 
was at x=0 and the only term left in Equation (2) was cos(wt). 
The reference signal for channel one was then reintroduced to 
the adder and its gain adjusted until the summed output signal 
was minimized while maintaining equality between the reference 
signal and the detector signal. The process was then repeated 


for channel two. 


C. THIRD ATTEMPT 

After the gains on the various circuit elements were 
adjusted, a second data set was collected. For this run, the 
mask was initially placed in a position so that when x=0, 
channel 1 was maximal. The results are listed in Appendix C 
and plotted in Figure 12. 

1. Problems 

7 The plot in Figure 12 shows that, although the 

amplitude of the periodic anomaly is reduced from its previous 


level, it is still present. It can be seen that the 
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Figure 12 Phase vs Position for the Second Data Set 


periodicity of the perturbation is still the same-.as the 
periodicity of the mask. The next cause of this periodicity 
to be investigated was a possible mechanical misalignment of 
the optical elements themselves. 
2. Investigation 

The potential for misalignment was examined by 
measuring the voltage output of the two detector channels. 
The mask was moved in increments of 0.1" and the channel 
voltages were read directly from the oscilloscope. | When 


plotted the data should be of the functional form 1+cos(x) for 
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channel one and i+sin(x) for evannel two, if the signal gains 
and mechanical alignment are perfect. The actual results are 
listed in Appendix D, while Figure 13 is a plot of the data. 
The expected curves are shown by the solid and dotted lines: 
the solid line indicates the curve 1+cos(x) while the dotted 
line indicates 1+sin(x) (both curves are modified to match the 
amplitude and offset of the data). It shows that, although 
channel one was approximately of the form 1+cos(x), channel 
two showed a significant departure from the curve 1+sin(x). 
This indicates that there was indeed a geometric problem with 


the setup of the optical elements. 





0 O01 02 O3 04 O5 O6 O7 O8 O98 1 = 21:1 
X Gin) 
—— 1+cos()} >-~ 1+sin(x) * Channel 1 4 Channel 2 


Figure 13 Detector Voltage vs Mask Position Prior to 
Realignment of Optical Elements 
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3. Solution 

The misalignment problems were easily reduced by 
realigning the optical elements as accurately as possible 
using calipers and a level. After the realignment was 
completed, the circuit amplifier gains were readjusted as 
described above in section IV.B.2. The test described in 
section IV.C.2 was run again and the results are plotted in 
Figure 14 (the data are listed in Appanale E). The solid line 
indicates A[1+cos(x)] and the dotted line Blitsin(x) ] were the 
constants A and B were added to match the amplitudes of the 
data. These results are much closer to expectations, thereby 
proving that precise mechanical alignment is mandatory for 


these types of experiments. 


D. FINAL EXPERIMENT 

Upon completion of the geometric realignment and 
adjustment of the amplifier gains, a final set of data was 
taken. A Mitutoyo Corp Code 500-351 digital caliper, accurate 
to +0.01lmm, was used to measure the displacement of the 
optical mask. Measurenents were made at mask displacement 
intervals of 0.5mm. The phase was once again measured with 
the lock-in amplifier. 

- The mask was displaced a total of 32mm and the caliper was 
checked several times to ensure that it stayed properly 
calibrated during the run. Also, the reproducability of the 


data was tested by moving the mask back to a position that had 
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Figure 14 Detector Voltage vs Mask Position After 
Realignment of the Optical Elements 


already been measured and comparing the phase difference with 
the previously measured value. ae was determined that the 
phase data was reproducible within +0.3°. The data are 
recorded in Appendix F and graphed in Figure 15. 

The periodic perturbation is still present, but it is not 
as pronounced as in preceding trials. It is believed that the 
ants is acceptable given the shortcomings of the experimental 


apparatus. 
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Figure 15 Graph of Phase vs Position for Final Data Set 


E. ANALYSIS 
1. Checking the Phase Relationship of the Driving Signals 
A possible source for the remaining periodic 
discrepancy was thought to be an extraneous phase shift 6 
between the two LED driving signals. These two signals were 
connected directly to the lock-in amplifier, and the phase 
difference between them was measured. The results of this 
test disclosed that the phase difference is in fact 


Ag=91.940.1° giving a value for 6 of 1.9°. 


35 











This extra phase 6 is of some concern. Because of its 
presence, the summation is not the desired summation as stated 
in Equation (1), but is rather 

sin(wt) sin(x) +cos (wt+8) cos (x). (3) 
By introducing the trig identity 
cos (wt+8) -cos (wt) cos (8) -sin(wt)sin(8), (4) 
to Equation (3), the second term can be rewritten as 
cos (®t) cos (x) cos (8) -sin(@t) sin (8) cos (x) (5) 
By using small angle approximations, sin(é)*6 and cos(6)#1, 
Equation (3) becomes, 
sin (®t) sin(x) +cos (®t) cos (x) -3sin(®t) cos (x). (6) 


To determine the overall effect of the extraneous 
phase shift 6, the measured signal is assumed to be cos(ét-A) 
where A is the actual measured phase difference. This can be 


rewritten as 
cos (wt-A)=-cos (wt) cos(A)+sin(wt)sin(A). (7) 


Comparing the right hand side of this equation to Equation (6) 


gives the following results, 


cos (A) =cos (x) 


(8) 


sin(A) «sin (x) -8cos (x) . 
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And finally, solving Equations (8) for A yields 


a sin (x) -8cos (x) - . (9) 


A-tan 
cos (x) 


If the residual electrical phase difference 6 was 
indeed the cause for the periodic vontineaeity in the final 
data set (Figure 15), then plotting A vs x should result ina 
Simulated phase curve similar te. the actual experimental 
result. Such a plot, with 6=1.9°, is shown in Figure 16. 
This graph does show a panisaie. anomaly, however it is 
symmetric about the ideal linear curve while the actual 
experimental data is not. Therefore, although the phase 
difference § must be a contributing factor in the periodic 
discrepancy in the measured data, it cannot be the only cause. 

2. Further Investigation 

Another possible cause is that the gain matching 
discussed in section IV.B.2. was not completely successful. 
This would mean that the sin(@t) and cos(®#t) terms were not 
fully subtracted out from Equation (2) and that the amplitudes 
of the summed signals in Equation (1) were not identical. If 
this is true, then the signal being measured, after the high 


pass filter, is 


Asin(®C) +Bcos (®t) +Csin(®@t) sin (x) +Dcos (®t) cos (x)+5sin(@t)cos(x) (10) 
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Figure 16 Simulated Phase Curve for 6=1.9° | < 
Which can be written as 
sin(@ 4 [A+Csin (x) +3cos (x) ]+cos (wt) [B+Dcos (x) ] (11) 


Assuming that this equation corresponds to cos(®@t-A), and 


comparing it to the right hand side of Equation (7) gives 


cos (A). =B+Dcos (x), (12) 


and 


sin(A) -A+Csin(x) +8cos (x). (13) 
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By dividing Equation (13) by Equation (12), it is possible to 
solve for A: 


A+ Csin (x) +8cos (x) (14) 


rane): _ B+Dcos (x) 


To determine the coefficients A, B, C, and D, it is. 
necessary to perform a least-squares fit on the data 
[Ref. 9]. To simplify matters, a least-squares fit on tan(A) 
vs tan(x) was attempted. Carrying out the differentiation on 
the function and plugging into the least-squares fit formula 


gives the following system of four equations in four unknowns: 


A+Bsin(x,)+8cos(x;) _ Fe . 
ancien (15) 


2 [A+Bsin (x,) +8cos (x,))sin(x,;) _ f£,sin (x,) ae (16) 


[C+Dcos (x;)? C+Dcos (x;) 


: 2 j 
>> A+Bsin(x,) +8cos (x;) D> £{A+Bsin(x,) +8cos (x;) <b (17) 


[C+Dcos (x,)}? [C+Dcos (x;) }* 
> [A+Bsin (x,) +8cos (x,) ’cos (x;) >> f {A+Bsin (x,) +6cos (x,) ]cos (x;) -o (18) 
7 [C+Dcos (x;) is 7 [C+Dcos (x;) FP 
It was determined in discussions between the author 
and his advisor that solving this system of equations would 
involve techniques that are beyond the scope of this thesis. 
In order to further investigate the effect of gain mismatch 
another attempt at performing a partial least-squares fit was 
made, using a subset of Equations (15)-(18). By setting the 
constants A and B in Equation (10) equal to zero, the gain 


mismatch where the sin(wt) and cos(wt) terms (Equation (2)) 
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are subtracted out is eliminated as a possible contribution to 
the phase anomaly. Equation (14) then reduces’ to 
tan(A)=atan(x)+8 where a=C/D and B=6/D. Solving the resulting 
systems of two equations and plotting A versus x gives the 
results in Figure 17. The displayed behavior of the curve 
_Goes not resemble the actual behavior of the measured data. 
Therefore it is concluded that any difference in amplitudes (C 
and D) between the sin(ot) sin(x) and the cos(wt)cos(x) terms 
in Equations (10) that may be present has negligible effect on 
the actual experiment. The observed anomaly is not produced 


by such amplitude mismatches. 





15 
X (mm) 


Figure 17 Simulated Phase Curve with Constants A=B=0 in 
Equation (10) 
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V. CONCLUSIONS /RECOMMENDATIONS 


A. CONCLUSIONS 


As stated in Chapter One, the goal of this thesis was to 
develop and to implement a proof-of-concept validation of the 
proposed electro-optical phase encoded position transducer. 
The plot in Figure 15 gives incontestable proof that the 
underlying concepts behind the technique are valid, although 
some residual problems remain. Further development is 
required before the technique can be incorporated into a 
working instrument for use in the infrared imaging system 


proposed by Davis. 


B. RECOMMENDATIONS 


The most difficult aspects of the project were the 
matching of the amplifier gains and the positioning of the 
optical elements. The mismatch of the amplifier gains was the 
most likely cause for the 1.9° phase error in the quadrature 
signals used to drive the LEDs. It probably also played a 
part in the nonlinearity of the final data. Possible error in 
the alignment of the optical elements may also have 


contributed to the observed nonlinearity. 
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It is recommended that, in future designs, higher 
precision electrical components be used in all circuits. This 
would eliminate the residual oe error as well as make 
matching the amplifier gains more stable and precise. 

A complete least-squares fit of the data using linear 
programming techniques on Equations (10) should be conducted 
to determine the significance of the amplifier gain problems 
in the final results, with corrective action taken as 
described above. 

The problem of alignment of the optical elements was a 
function of the equipment available to the author. Although 
great -care will need to be taken in future designs to ensure 
the alignment is accurate, this should not be a major problem 
when using more precise opto-mechanical micropositioning 


equipment. 
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APPENDIX A 





/*. Program THESPLOT 
‘ Written by Lt Gary R Parriott, USN 


This program will draw an optical mask consisting of a 
sine wave and a cosine wave such that the sine wave is two 
inches above the cosine. The areas under the curve will be 
white, while the background areas (above the curves) will be 
black. The mask will be surrounded by a black window. ¥*/ 


#include <stdio.h> 
#include <math.h> 


int main(void) 


{ 
/* Variable Declaration*/ 


float A,x,x0,y,xscl,yscl; 
double PI=3.14159; 

int x1,y0,y1,y2,y3,y4,period; 
i FILE *out; 


/* Open file named mask.dat in default directory and 
return error if not allowed */ 


if ( (out=fopen("mask.dat", "wt") )==NULL) printf ("No 
File/n") ; 


/* Enter the amplitude and number of periods required. */ 


printf("Enter the number of periods desired.\n"); 
scanf("%d", &period); 

printf("Enter the amplitude of the Sine waves in 
centimeters. \n"); 

scant ("%Zf", &A); 


/* Determine the values of the variables not designated by 
the user */ 


xsc1=4400/(2*period); /* Scales the X coordinate */ 
yscl=A*400; /* Scales the Y coordinate */ 
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y1i=6184; /* Sets the upper edge of the window */ 
y2=5000-yscl-10; /* Sets the lower edge of the Sine */ 
y3=2968-yscl-10; /* Sets the lower edge of the Cosine */ 
y4=1784; /* Sets the lower edge of the window */ 


/* Initialize plotter, set up initial parameters */ 
fprintf (out, "%s","IN;PTO.7;SP1;GM12752;") ; 7 f 
fprintf (out,"%s","PU3000,5000;") ; 

fprintf(out,"%s","PMO;PA;PD;") ; 

/* Determine coordinates of Sine wave for polygon */ 


for (x=0;x<=4400;x=x+25) 


¢ 
x0=x*PI/xscl; 
y=sin(x0) ;x1= (int) x+3000; yo=(int) (yscl*y) +5000; 
fprintf (out, "Sstdtstdzs","PD",x1,",",yO,";") ; 

} 


/* Define, fill and plot rest of polygon */ 

fprintf (out, "%sta%std%s","PD7400,",y1,";PD3000,",y1,"3") ; 
fprintf (out,"%s","PM2;FP;EP;") ; 

/* Reset parameters for new polygon */ 

fprintf (out, "%s","PU3000,2968;") ; 
fprintf(out,"%s","PMO;PA;PD;") ; 

/* Determine coordinates of Cosine wave for polygon */ 


for (x=0;x<=4400;x=x+25) 


{ 
xO=x*PI/xscl; 
y=cos (x0) ;x1= (int) x+3000;y0= (int) (yscl*y) +2968; 
fprintf (out, "Sstdtstdts","PD",x1,",",yO,";") ; 

} 


/* Define, fill and plot rest of polygon */ 


fprintf£ (out, "Sstd%s%d%s", "PD7400,",y2,";PD3000,",y2,";") ; 7 
fprintf(out,"%s","PM2;FP;EP;") ; 
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/* Define and fill rectangles to finish plot */ 

fprintf (out, "Ss%d%s%ds", "PU2500,",y4,";RA3000,",y1,"3");7 
fprintf (out, "%s%d%s%d%s", "PU3000,",y4,";RA7400,",y3,"7");7 
fprintf (out, "%s%d%s%da%s", "PU7400,",y4,";RA7800,",y1,"7");7 
/* End plotting routine */ 


fprintf(out,"%s","SP0;PG;"); 


/* Close file and end program*/ 


fclose (out) ; 
return; 
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APPENDIX B 


Phase vs Position for the Initial Data Set 
Plotted in Figure 11 
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APPENDIX C 


Phase vs Position for the Second Data Set 
Plotted in Figure 12 
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APPENDIX D 


Detector Voltage vs Mask Position 
Prior to Realignment of Optical Elements $ 
Plotted in Figure 13 
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APPENDIX E 


Detector Voltage vs Mask Position 
After Realignment of Optical Elements 
\ Plotted in Figure 14 
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APPENDIX F 


Phase vs Position 
Final Data Set 
Plotted in Figure 15 t 
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